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Abstract 

Using a sample of 102 million T(15) events collected with the Belle detector, we report on 
the first search for charge-parity-even charmonium and charmonium-like states in T(1S) radia- 
tive decays. No significant Xcj or n c signal is observed and 90% C.L. limits on B(T(1S) — > 
7Xco) < 6.5 x lO" 4 , B(T(1S) -> 7 Xci) < 2.3 x 10~ 5 , B(T(1S) -> 1X c2) < 7.6 x KT 6 , and 
B(T(1S) — > 777c) < 5.7 x 10~ 5 are obtained. The product branching fraction limits B(T(1S) — > 
7 X(3872))S(X(3872) -»• 7r+7r"J/^) < 1.6 x 10~ 6 , B(T(1S) -> 7 X(3872))fi(X(3872) -»• 
7r + 7r-7r J/V>) < 2.8 x 10~ 6 , B(T(15) -»• 7 X(3915))£(X(3915) -> uoJ/ip) < 3.0 x 10~ 6 , and 
B(T(1S) -»• 7y(4140))/3(y(4140) -»• cpJ/ip) < 2.2 x 10~ 6 are obtained at the 90% C.L. Fur- 
thermore, no evidence is found for excited charmonium states below 4.8 GeV/c 2 . 

PACS numbers: 14.40.Nd, 14.20.Lq, 13.25.Gv 
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There is renewed interest in charmonium spectroscopy after the operation of the two 
5-factories. In addition to many conventional charmonium states, a number of states with 
unusual properties have been discovered, which may include states beyond the quark-model, 
such as quark-gluon hybrids, meson molecules, multi-quark states, and so on hj-m. States 
with J PC = 1 can be studied using initial state radiation (ISR) in the large T(4S) data 
samples. For the study of charge-parity-even charmonium states, radiative decays of the T 
states below open-bottom threshold are used. 

The production rates of the lowest lying P-wave spin-triplet (xcj, J=0, 1, or 2) and S- 
wave spin-singlet (i] c ) states in T(1S) radiative decays are calculated in Ref. |8j, where the 
former is at the part per million level, and the latter is about 5 x 10 -5 . There are no calcula- 
tions for radiative decays involving excited charmonium states, let alone for charmonium-like 
states, such as the X(3872) Q, the X(3915) @, and the F(4140) Q. 

In this paper, we report on a search for the Xcj, Vc X(3872), X(3915), and F(4140) 
states in Y(IS') radiative decays. The Xcj states are reconstructed via their El transition to 
the J/ip. The 1] C is reconstructed in the KgK + Tc~ + c.c, tc + it~K + K~, 2(K + K~), 2(7r + 7r~), 
and 3(7r + 7r~) final states. To search for the X(3872) and X(3915), we use the tt + ti~ J/ip 
and n + n~ n° J / ip final states, while we reconstruct the F(4140) in the 4>J/ip mode. This 
analysis is based on a 5.7 fb _1 data sample collected at the T(IS') (102 million Y(IS') 
events) and a 1.8 fb _1 data sample collected at \fs = 9.43 GeV (continuum data) with the 
Belle detector [Io[ operating at the KEKB asymmetric-energy e + e~ collider (Tlj . 

For each charged track, the impact parameters perpendicular to and along the beam 
direction with respect to the interaction point are required to be less than 0.5 cm and 4 cm, 
respectively, and the transverse momentum must exceed 0.1 GeV/c in the laboratory frame. 
For each charged track, information from different detector subsystems is combined to form 



a likelihood for each particle species M2|. Tracks with TZ K = r L l r > 0.6 are iden- 



tified as kaons, tracks with IZk < 0.4 are identified as pions. With these selections, the 
kaon (pion) identification efficiency is about 90% (96%), while 5% (6%) of kaons (pions) are 
misidentified as pions (kaons). For electron identification, the likelihood ratio is defined as 
T^e = r ) where C e and C x are the likelihoods for electron and non-electron hypotheses, 
respectively. These are determined using the ratio of the energy deposited in the electromag- 
netic calorimeter (ECL) to the momentum measured in the silicon vertex detector (SVD) 
and central drift chamber (CDC), the shower shape in the ECL, the matching between the 
position of the charged track trajectory and the cluster position in the ECL, hit informa- 
tion from the aerogel threshold Cherenkov counters (ACC), and dE/dx information in the 



CDC LL3J. For muon identification, the likelihood ratio is defined as TZ a = ~ , r M , r , where 



£ M , C n and Ck are the likelihoods for muon, pion and kaon hypotheses, respectively. These 
are based on track matching quality and penetration depth of associated hits in the iron 
flux-return (KLM) [ll]. 

We reconstruct J/ijj mesons from e + e~ or candidates. In order to reduce the effect 

of bremsstrahlung or final-state radiation, photons detected in the ECL within 0.05 radians 
of the original e + or e~ direction are included in the calculation of the e + e~(7) invariant 
mass. For electrons from J/ip — > e + e~, one track should have lZ e > 0.95 and the other 
7l e > 0.05; for muons from J/ip — > at least one track should have TZ^ > 0.95 (in the 

Xcj analysis, the other track should have associated hits in the KLM detector that agree 
with the extrapolated trajectory of a charged track provided by the drift chamber). The 
lepton identification efficiency is about 90% for J/ip — > e + e~ and 87% for J/ijj — > ji + . 
In order to improve the J/ip momentum resolution, a mass fit to the reconstructed J/ijj 
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candidates is then performed for all the channels with J/ip signals. 

A neutral cluster is considered to be a photon candidate if its ECL shower does not match 
the extrapolation of any charged track and the energy deposition is greater than 40 MeV. 
The photon candidate with the maximum energy in the e + e~ center-of-mass (CM.) frame 
is taken to be the T(1S) radiative decay photon, and its energy is required to be greater 
than 3.5 GeV, which corresponds to a 4.8 GeV/c 2 mass particle produced in T(1S) radiative 
decays. 

To study the 7Xcj mode, we reconstruct Xcj via its decay into jJ/ip. The deposited energy 
of Xcj's photon is required to be greater than 150 MeV, and the total number of photons 
in the event is required to be exactly two, in order to suppress multi-photon backgrounds. 
The higher energy photon is denoted as 7^ and the lower energy one is denoted as 7;. The 
angle between the two photons should be larger than 18° to remove the background from 
split-off fake photons. To remove the ISR background e + e~ — > 7isr'0(2S) — >■ 7isr7XcJ, 
where a photon is missed, we require the square of the missing mass of 7/ and lepton-pair 
to be within -0.5 GeV 2 /c 4 and 0.5 GeV 2 /c 4 since this background has at least two missing 
photons (the 7isr photon(s) and one photon from the ip(2S) decay) and the missing mass 
tends to be large. Bhabha and dimuon background events with final-state radiative photons 
are further suppressed by removing events where a photon is detected within a 18° cone 
around each charged track direction. 

A clear J/ip signal is observed in the mode, while no significant J/ip signal is 

observed in the e + e~ mode due to residual radiative Bhabha background. The J/ip signal 
region is defined as \mg+£- —mj/^\ < 30 MeV/c 2 (~ 2.5a), and the J/ip mass sidebands are 
defined as 2.959 GeV/c 2 < m e+e - < 3.019 GeV/c 2 or 3.175 GeV/c 2 < m £+£ - < 3.235 GeV/c 2 ; 
where the latter is twice as wide as the signal region. 

Figured] shows the 'jiJ/ip invariant mass distribution after the above selections are applied 
to the T(15) data sample for the combined e + e~ and /i + /U~ modes, together with the 
background estimated from the normalized J / ip mass sidebands. Apart from possible weak 
Xco and Xci signals, the J / ip sideband events represent well the signal region, indicating that 
the production of any of the \ c j states is not significant. There are no structures at higher 
masses, where we would expect excited Xcj states. 
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FIG. 1: The ~fiJ/ip invariant mass distribution in the T(1S) data sample. Hints of Xco and 
Xci signals are seen although no obvious Xc2 signal is observed. The solid curve is the best fit, 
the dashed curve is the background, and the shaded histogram is from the normalized J/ip mass 
sidebands. 



A fit to the signal region is performed with Breit-Wigner (BW) functions for the reso- 
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nances convolved with Gaussian resolution functions and a second-order polynomial back- 
ground term. This fit yields 5.9i|i, 8.5t|;f, and 0.6^4 events for the Xco, Xa, and Xc2, 
respectively. Here the width of the Gaussian resolution function is fixed as 7.0 MeV/c 2 , its 
MC-determined value. Bayesian upper limits on the number of events at the 90% CL. by 
integrating the likelihood distribution (as a function of the yield) are found to be 11.5, 13.8, 
and 2.4 for the Xco, Xci, an d Xc2, respectively. 

To study the 77/ c mode, we reconstruct the r] c mass from the invariant masses of 
KgK + ir~ + c.c, n + n-K + K-, 2(K + K~), 2(7T + 7T~), and 3(n + n-). Well-measured charged 
tracks are selected and the numbers of charged tracks are six for the 3(7r + 7r~) final state 
and four for the other final states. All the charged tracks are required to be identified as 
kaons or pions. The recoil mass-squared of the charged particles in each 7] c decay mode is 
required to be within —1 GeV 2 /c 4 and 1 GeV 2 /c 4 . For Kg candidates decaying into 7r + 7r~ 
in the K s K + tt~ + c.c. mode, we require that the invariant mass of the 7r + 7r~ pair lie within 
30 MeV/c 2 of the K° s nominal mass and that the Kg candidate must have a displaced vertex 
and flight direction consistent with a Kg originating from the IP; the same selection method 
is used in Ref. [IBj]. There are events with leptons misidentified as pions in the 7r + 7r _ K + K~ 
and 2(7r + 7r~) modes, and they are removed by requiring lZ e < 0.9 and IZ^ < 0.9 for the pion 
candidates. 

Figure |2] shows the combined mass distribution for the five r\ c decay modes after the 
selection described above. The peak in hadronic mass at the J/ip mass, as seen in Fig. [2j 
can be attributed to the ISR process, e + e~ —> ■JisrJ/4', while the accumulation of events 
within the r) c mass region is small. The shaded histogram in Fig. [2] is the same distribution 
for the continuum data, normalized according to the ratio of the luminosities on and off the 
T(IS') peak. From Fig. 121 we can see that the J ftp signal in T(IS') data is well reproduced by 
the normalized continuum data, demonstrating its ISR origin. It is also evident that T(1S) 
radiative decays to light hadrons are substantial, as indicated by the difference between the 
number of non-resonant events in the two data sets. 
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FIG. 2: The mass distribution for a sum of the five rj c decay modes. The solid line is a sum of 
the corresponding functions obtained from a simultaneous fit to all the rj c decay modes, and the 
dashed line is a sum of the background functions from the fit. The shaded histogram is a sum of 
the normalized continuum events, where the J/ip signal is produced via ISR. 



We perform a simultaneous fit to all the r\ c decay modes, where the r] c mass and width 
are taken from the PDG [l6j], and the ratio of the yields in all the channels is fixed to Bj€i, 



where each £>j is the rj c decay branching fraction for the z-th mode reported by the PDG [16 



and e, is the MC-determined efficiency for this mode. In the fit, we take a BW function 
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convolved with a Gaussian resolution function (its resolution is fixed to 7.9 MeV/c 2 from 
MC simulation) as the r] c signal shape, another Gaussian function as the J/ip signal shape, 
and a second-order polynomial as the background shape. The fitted results are shown in 
Fig. [21 where the solid line is the sum of the best fit functions in the simultaneous fit, and 
the dashed line is the sum of the background functions. The fit yields 46 ± 22 r] c signal 
events, with a statistical significance of 2.2cr. An upper limit on the number of the r\ c signal 
events is estimated to be 72 at the 90% C.L. From the fit, we obtain 89 ± 20 and 54 ± 16 
J/ijj signal events in T(IS') and normalized continuum data samples, respectively, with a 
mass of 3099.9 ± 2.1 MeV/c 2 , which is consistent with PDG value. 

The selection criteria for T(IS') — > 7X(3872h X(3872) — > 7i + 7i^J/ip are similar to those 
used for ISR tc + tc~J/i/j events in T(AS) data Q. We require that one J J if) candidate be 
reconstructed, that two well identified 7r's have an invariant mass greater than 0.35 GeV/c 2 , 
and that the recoil mass-squared of the TT + n~J/tp be between —1 GeV 2 /c 4 and 1 GeV 2 /c 4 . 
To suppress ISR 7i + n~J/tp background, we require that the polar angle of the radiative 
photon satisfy | cos 6*| < 0.9 in the e + e~ CM. system. Except for a few remaining ISR 
produced '0(25') signal events, only a small number of events appear above the ip(2S) peak 
in the n+n^J/ip invariant mass distribution, as shown in Fig. [3]^a). Within the X(3872) 
signal region, there is one event with a mass of 3.870 GeV/c 2 . However, there are no events 
in the J/ip mass sidebands from 3.6 to 4.8 GeV/c 2 . We estimate the statistical significance 
of the X(3872) signal to be 2.3cx if the background distribution is flat above 3.7 GeV/c 2 . 
Assuming that the number of signal events follows a Poisson distribution with a uniform 
prior probability density function and there is no background, the upper limit on the number 



of the X(3872) signal events is 3.9 [16 
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FIG. 3: (a) Distribution of the ir + tt~ J j 'ip invariant mass for Y(1S) — > jTr + n~ J / tp candidates, 
(b) Distribution of the tt + tt~tt° J/tp invariant mass for T(IS') — > 77r + 7r~7r° J/ip candidates. Points 
with error bars are data, open histograms are the MC expectation for the X(3872) signal (not 
normalized). The peak at 3.686 GeV/c 2 in (a) is due to ip(2S) production via ISR. 



We validate our analysis by measuring the ip(2S) ISR production cross section as observed 
in the 7r + 7r~J/?/> mode. By relaxing the photon polar angle requirement, we observe 383 
ip(2S) signal events and a cross section of e + e~ — > 7 ISR , 0(2S') is measured to be (20.2 ± 
1.1 (stat.)) pb, in agreement with a theoretical calculation of 18.5 pb using PDG [l6[ values 
for the ip{2S) resonance parameters as input. 

To study the / ~f7T + 7i^7i J/ifj mode, we require the invariant mass of a pair of photons to be 
within 10 MeV/c 2 around the nominal 7r° mass (the mass resolution is about 4 MeV/c 2 ) to 
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select 7T° candidates. The other event selection criteria are similar to those in the X(3872) — > 
7T + 7T~J/ip mode, except that we do not require the 7t + tt~ invariant mass to be greater than 
0.35 GeV/c 2 , which is used to remove the 7 conversion background events in the 7t + tt~J /ip 
mode. 

Figure El^b) shows the 7r + 7r~7r° J/ip invariant mass distribution, where the open histogram 
is the MC expectation for the X(3872) signal shape plotted with an arbitrary normal- 
ization. We observe two events in the 7r + 7r~7r° J/-?/> mass spectrum between 3.6 GeV/c 2 
and 4.8 GeV/c 2 in the T(IS') data. For these two events, the 7r + 7r~7r° J/ijj masses are 
3.67 GeV/c 2 and 4.23 GeV/c 2 , and the corresponding tt + 7t~7t° masses are 0.54 GeV/c 2 
and 1.04 GeV/c 2 , respectively. The event at 3.67 GeV/c 2 , is likely to be from e + e _ — > 
'jiSRTjJ/ip — > 7isR7r + 7r _ 7r°^ + ^~, since 0.9 events are expected from MC simulation. No event 
is observed within the X(3872) or X(3915) mass region. An upper limit on the number of 
X(3872) or X(3915) signal events is 2.3 at the 90% C.L. 0. 

We also search for the F(4140) in its decays into 4>J/i/j, with cf) — > K + K~~ and 
J/ip — > £ + £~. The selection criteria are very similar to the analysis of X(3872) — > 7r + 7r~ J/ip 
above. Here two kaons are required to be positively identified and one J/ip candidate is 
reconstructed. No clear J/ip or signal can be seen after the initial event selection. We 
define the J/ip signal region as \mi+g- — mj/$\ < 30 MeV/c 2 , and the signal region as 
1.01 GeV/c 2 < m K + K - < 1.03 GeV/c 2 , according to MC simulation. After applying all of 
the above event selection criteria, there are no candidate events in the <j)J/ip invariant mass 
region between 4 GeV/c 2 and 4.8 GeV/c 2 . An upper limit on the number of F(4140) signal 
events is 2.3 at the 90% C.L. fil. 



There are several sources of systematic error in determining limits on the branching 
fractions. A particle identification efficiency uncertainty between 2.4%-3.7% is assigned 
depending on the final state particles. An uncertainty in the tracking efficiency for tracks 
with angles and momenta characteristic of signal events is about 1% per track, and is 
additive. Photon reconstruction contributes an additional 2% per photon. Errors on the 
branching fractions of the intermediate states are taken from the PDG [HI]. F° r the 7] c decays, 
the biggest difference in the efficiency by using a phase space distribution and including 
possible intermediate resonance states is 2.1%. The difference in overall efficiency for a 
flat radiative photon angular distribution and a 1 ± cos 2 6 distribution is less than 3.0%. 
Therefore, we quote an additional error of 5% for all the states studied due to limited 
knowledge of their decay dynamics. According to MC simulation, the trigger efficiency 
is rather high, with an uncertainty that is smaller than 1%. The uncertainty due to the 
missing mass squared requirement is 1.0% for the channels with only one photon and 4.7% 
for channels with more than one photon. Uncertainties on the Xcj arid r] c signal event yields 
are estimated to be 1.6% and 15%, respectively, by changing the order of the background 
polynomial, the range of the fit, and the values of the masses and widths of the resonances. 
In the T(IS') — > 'JXcj mode, the uncertainty that is associated with the requirement on the 
number of photons is 2% after applying a correction factor of 0.96 to the MC efficiency, 
which is determined from a study of a very pure T(IS') — > fi + fi~ event sample. In the 
i] c — > KgK + ir~ + c.c. mode, the uncertainty in the K$ efficiency is determined by comparing 
yields for a sample of high momentum K$ — > 7i + tt~ decays before and after applying the 
Ks candidate selection criteria; the efficiency difference between data and MC simulation 
is less than 4.9% [13]. Finally, the uncertainty on the total number of T(IS') events is 
2.2%. Assuming that all of these systematic error sources are independent, we add them 
in quadrature to obtain total systematic errors as shown in Tab. [B In order to calculate 
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conservative upper limits on these branching fractions, the efficiencies have been lowered by 
a factor of 1 — cr sys . 

In summary, Table [J lists the final results for the upper limits on the branching fractions 
of all the states studied, together with the upper limits on the numbers of signal events and 
their detection efficiencies. The results obtained on the \cj and rj c production rates are not 
in contradiction with the calculations in Ref. jgj]. No X(3872), X(3915), or F(4140) signals 
are observed, and the production rates of the n + 7i~ J/ip, n + n~ n J / tp , uJ/tf), or <pj '/if) modes 
are found to be less than a few times 1CT 6 at the 90% CL. Furthermore, we find no evidence 
for excited charmonium states below 4.8 GeV/c 2 . 

TABLE I: Summary of the limits on T(IS') radiative decays to charmonium and charmonium-like 
states R. is the upper limit on the number of signal events, e is the efficiency, <r sys is the total 
systematic error and B(T(IS) — > 7i?) UP (Br) is the upper limit at the 90% C.L. on the decay 
branching fraction in the charmonium state case, and on the product branching fraction in the 
charmonium-like state case. 
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